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ABSTRACT 



We report optical imaging, optical and near-infrared polarimetry, and Spitzer 
mid-infrared spectroscopy of comet C/2007 N3 (Lulin). Polarimetric observa- 
tions were obtained in R (0.676 /im) at phase angles from 0.44° to 21° with 
simultaneous observations in H (1.65 /im) at 4.0°, exploring the negative branch 
in polarization. Comet C/2007 N3 (Lulin) shows a typical negative polarization 
in the optical as well as a similar negative branch at near-infrared wavelengths. 
The 10 /im silicate feature is only weakly in emission and according to our ther- 
mal models, is consistent with emission from a mixture of silicate and carbon 
material. We argue that large, low-porosity (akin to Ballistic Particle Cluster 
Aggregates), rather absorbing aggregate dust particles best explain both the po- 
larimetric and the mid-infrared spectral energy distribution. 

Subject headings: comets: general — comets: individual (C/2007 N3 (Lulin)) 
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1. INTRODUCTION 

Comet dust provides us with insight into two fundamental aspects of our solar system: 
(a) the composition of both the solar system's initial ingredients and formation products; 
and, (b) the spatial and temporal heterogeneity of the proto-planetary disk. Comet dust 
comae are optically thin, providing a unique laboratory for the study of preserved dust 
grains that observations of Kuiper Belt objects and primitive asteroids cannot provide. The 
compositions of both comet dynamical families, ecliptic comets (ECs) and nearly-isotropic 
comets (NICs), have been preserved in their interiors during their residences in the Scattered 
Disk and Oort Cloud. 



Knowledge of cometary materials enables better understanding of planetary and small 
body formation as well as evolutionary processes in the protoplanetary disk. The NASA 
space mission Stardust was designed to bring back samples of comet dust from 8 IP/ Wild 
2. The terrestrial laboratory analysis of grain samples captured in the aerogel have helped 
us to understand the properties of co met refractor y materials, revealing, among others, 



high-temperature crystalline materials ( IFlynn 



20081 ). In the preliminary analysis of Stardust 



samples, ~ 65% of the aerogel tracks were 'carrot-shaped,' indicative of large solid grains 



that upon inspection inclu ded micron- to tens o 



silicate) and FeS crystals (IBurchell et al. 



2008 



micron-sized 



Zolensky et al. 



orster ite (Mg-rich crystalline 
2008k A fraction, ~ 35%, 



of the Stardust aerogel tracks were bulbous and attributed to porous aggregates which 



suffered significant heati ng upon impact with the aerogel ( iRoskosz. Leroux. fc Watson 



2008 



Leroux et al. 



20081 ). However, analysis of microcrater impacts in the foil suggest that 



85% are clustered with irregular, non -circular outlines likely originating from composite, 



i.e., aggregate grains (IBorg et al. 



20071 ). In contrast to the Stardust samples that were 
dominated by weakly bound aggregates, which likely originated in the inner protoplanetary 
disk, the interplanetary dust particles (IDPs) captured in Earth's stratosphere from comet 
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26P/Grigg-Skjellerup were all highly porous (e.g., iNakamura- Messenger et al.ll2008l ) and 
were not dominated by large single mineral grains. 

Grain thermal modeling of mid-infrared (MIR) remote sensing spectral energy 
distribution (SEDs) of comets constrains various coma grain properties including their: 
temperatures, size distributions (n(a)da ), mass, dust production rate, and porosity (e.g., 



Harker et al 



2002 



Wooden et al 



19991 ). These grain properties must then be reconciled 



with those derived from analysis of the optical and near-infrared (NIR) polarimetry to 
produce a self-consistent description of the dust. 



'he status of remote sensing of cometary dust was reviewed by 



Kolokolova et al. 



( 120041 ). who concluded that aggregate models are likely best at exp laining photometric, 



polarimetric, and thermal properties of cometary dust. This model ( jKolokolova et al. 



20071 ) 



not only explained the phase dependence of brightness and polarization f or comets, but also 



explai ned t he existence of two classes of com e ts based on polari metry 



jl986|) and 



Levasseur-Regourd et al. 



J199J). 



Lisse et al. 



(|2002j) and 



by 



Dobrovo 



Sitko et al. 



sky et al. 



( 120041 ) 



argue that there are two classes of comets based on the strength of the silicate feature 
in emission and the 10 fim SED temperature excess over an equivalent blackbody at the 
observed heliocentric distance (r^) of the comet. Taxomomically, comets with strong silicate 
feature emission comprise the same population of comets that show high polarization 
(coma dust population dominated by small or very porous particles), while those with 
weak silicate feature emission are the same population that exhibit low polarization (coma 



dust population dominated by large compact 



with aggregate dust models (jKolokolova et al. 



particles). This correlation is consistent 



20071 ). Comparison of thermal model dust 



parameters used to interpret ground-based remote sensing observations of the SED arising 
from coma grains with the dust from the Stardust sample return mission has helped to 
improve the aggregate model, making it not only qualitatively but also quantitatively 
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compatible with the experimental data (IKolokolova &: Kimurall2010f ). 



The presence of the so-called negative branch in the phase depe ndence of the polarization 



is an important o 



Kolokolova et al. 



jserv ational characteristic of comet polarization (IMishcnenko et al. 



2010; 



2004 ). Unlike Rayleigh particles, which polarize light perpendicular 
(positive Stokes Q) to the plane of scattering at all phase angles, for comets and asteroids 
the p olarization position angle flips by 90° (negative Stokes Q) at low phase angles, < 20 ° 



e.g., 



Kelley et al. 



20041 ) . Although negative polarization can be rep roduced by mode 



comet dust as an ensem ble of polydisperse spherical particles (e.g., iMukai et al.lll987l ). 



mg 



Kolokolova et al. 



(120041 ) demonstrate that all other polarimetric, photometric, and thermal 
IR characteristics of comet dust require that comet partciles be aggregates of submicron 
monomers. This aggregate "s tructure" wa s confirmed by the in situ coma particle sampling 



of the Stardust mission (e.g., 



Flynn 



20081 ) that showed that the majority (~ 65%) of the 
returned sample particles had an aggregate nature. Reproduction of the negative branch 
can be achi eved in the aggregate m odel cases utilizing large (thousands of monomers), 



aggregates (IKolokolova et al 



transparent silicates ( 



2007!) or 



Zubko et al. 



2009 



alternatively i n cases where the aggregates are 



made of transparent silicates is hard to reconci 
nuclei caused by a sig nificant carbon content ( 



Kimurall200ll ). However, the presence of aggregates 
e with the low geom etric albedo of comet 



aggregate models (e.g., 



Lasue et al. 



2009 



Campins et 



Shen et al 



2009 . 



al 



20061 ). Other ballistic 



2008) developed to assess the 



effect of effective aggregate size, porosity, and composition in the polarimetric behavior of 
dust in astrophysical environments need to be carefully scrutinized for their applicability 
in describing dust in comet coma. This assessment, which we discuss in detail in §3.41 is 
critical for interpreting the negative polarization branch observed at small phase angles. 

There are very few observations of comet polarization in the NIR at phase angles below 
20°, and existing observations present a puzzle. The unusual comet C/1995 01 (Hale-Bopp), 
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for example, showed no negative branch at NIR wavelengths (I Jones fc Gehrzll2000l ). unlike 
the behavior of all comets at optical wavelengths. Knowledge of the wavelength dependence 
of the polarization at long wavelengths (A > 1.3 /xm) through the negative branch as well 
as NIR observations at a wide range of phase angles should shed more light on the nature 
of the aggregates that make up the majority of comet dust. These aggregates also need to 
reproduce the observed emission in the thermal infrared, especially in the 10 /xm spectral 
region where various minerals produce broad resonances and narrow features. 



The wavelength dependence of the pol arization has on l y rece ntly 



Deen extended out 



beyond 1 /xm. By observing out to 2.2 /im, 



Jones fc Gehrzl f l2000h and 



Kelley et al 



(120041 ) 



were able to show that in general, the polarization increases from the optical out to the 
NIR for most comets (red polarization color). This effect is in contrast to that observed 
in asteroids, where the polarization decreases with wavelength (blue polarization color; 



Masiero fc Celling 



2009 



, and references therein). Most NIR polarimetry observations of 
comets are made independent of the visual observations, requiring extrapolation between 
different phase angles in order to directly compare fractional polarization at visual and NIR 
wavelengths. This lack of simultaneity complicates determination of the polarization color. 

Comet C/2007 N3 (Lulin) provided an excellent opportunity to make simultaneous 
visual and NIR polarimetric observations of a comet at phase angles below 20°. Comet 
C/2007 N3 (Lulin) is a NIC (Oort cloud) comet that passed less than 1° from zero phase 
angle at opposition (see Fig. [1]). We obtained polarimetric observations in the optical red 
at phase angles from 0.44° to 21° with simultaneous observations in the NIR at 4.0°. When 
combined with Spitzer MIR spectroscopy, we can compare dust models that match the 
observed scattering and the observed emission properties of the dust in the coma. 
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OBSERVATIONS 



2.1. Optical Imaging 

Comet C/2007 N3 (Lulin) was observed on 2009 March 03 UT with the 2.3-m Bok 
Telescope at Kitt Peak National Observatory at heliocentric distance, = 1.45 AU and a 
geocentric distance, A = 0.49 AU. The comet was at a phase angle a = 17.77°. 



20041 ). a prime 



The images were obtained with the 90prime Camera ([Williams et al. 
focus imager built for the Bok Telescope. At the time of observation, the 90prime camera 
utilized a thinned back-illuminated CCD detector with 4064 x 4064 pixels with a pixel 
size of 15.0 fim. At prime focus the camera pixel scale is 0.45" which yields a field of view 
of 30.5x30.5 square-arcmin. The instrument was equipped with Cousins/Bessel system 
broadband V and / filters. 

Multiple exposures (3 images of 120 seconds in I band and 4 images in V band; 1 of 
60 seconds and 3 of 120 seconds) were obtained of the nucleus and coma of the comet with 
the telescope tracking at the non-sidereal rate corresponding to the predicted motion of 
the comet provided by JPL Horizons j in an airmass range of 1.08 to 1.11. All images were 
corrected for overscan, bias and flat-fielding with standard IRAtj^ routines. Due to the 
thinned nature of the 90prime detector, / band data must be defringed. A master fringe 
file was created using the median of 41 data frames from this night of observation and 
subtracted from the / band data. The data was photometrically calibrated using Landolt 



standard fields PG0918, SA102, SA104, PG1323 flLandoltl 



1992l ) and the average nightly 



x http: / / ssd.jpl.nasa.gov/horizons.cgi 

2 IRAF is distributed by the National Optical Astronomy Observatory, which is operated 
by the Association of Universities for Research in Astronomy (AURA) under cooperative 
agreement with the National Science Foundation. 
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seeing was l'/8. One 120 second exposure in each band is shown in Fig. [2j 



2.2. Polarimetry 



The NIR obs ervations were made in the if(1.65 /im) band using NSFCAM2 



(jShure et al. 



1994 ) on the IRTF 3-m telescope at a plate scale of 0'.'04 per pixel, resulting 



in a field of view of 80" square. In polarimetry mode, NSFCAM2 utilizes a rotating 
half-wave plate at the entrance window of the camera and a cold wire grid polarizer in the 



second filter wheel. Details of the NSFCAM2 
reduction procedure are given in 



Jones fc Gehrzl (12000 ) and 



olari mete r) observing; 



polarimetric observations were made using OptiPol (Uones et al 



t echniq ue and data 



Kel lev et al. 



(12004 ) . The visual 



20081 ) at the University of 



Minnesota Mt. Lemmon Observing Facility (MLOF) in a narrow band R filter centered 
at A G = 0.676 /im, with a bandwidth AA = 0.04 /im. This filter was chosen to avoid 
contamination from unpolarized gas emission as much as possible. 

The results of our polarimetry, including the observing log, observation time, and 
phase angle are given in Table [TJ Simultaneous polarimetry at R and H was obtained on 
2009 January 27 UT. 



2.3. Spitzer IRS 



Comet C/2007 N3 (Lulin) was observed by Spitzer on 2008 October 04.3 UT 
[r^ = 1.90 AU, a Spitzer-co met distance of 1.6 73 AU, and a phase angle of 32.5°) with the 



Infrared Spectrograph (IRS; 



Houck et al 



20041 ) as part of a larger Cycle 5 study assessing 



the water production and volatile production rates of comets (program identification 
number [PID] 50335; PI: D. E. Harker). The astronomical observation request (AOR) key 
for the dataset is 2587584, and the basic calibrated data (BCD) products processed with 
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IRS reduction pipeline S18.60. The AOR for the short-wavelength, low-resolution SL1 
(7.4 — 14.5 /im) data discussed here executed a 7 x 3 spectral map (performed with no 
peak-up) yielding 21 spectra (6 sec x 2 cycles) with 1.87" x 10.0" steps (perpendicular x 
parallel to the long slit dimension). Background (shadow) observations were taken 34 hours 
later at the same celestial coordinates as the target spectra (AOR key 25988864), allowing 
the comet to move out of the sp ectral map field-q f -view (FOV). Further analysis of the SL1 



and LH spectra are discussed in 



Woodward et al. 



(2011) 



The spectra were reduced as follows. The shadow observations were subtracted from 
the on-source obser vations, and the r esult was assembled into a data cube using the 



CUBISM software flSmith et al. 



20071 ). with bad pixels masked and all extended source 
calibrations applied. The spectrum presented in this paper is extracted from an aperture 
9"25 x 9"25 in size, centered on the peak surface brightness of the comet. 



3. DISCUSSION 



3.1. Optical Imaging 

The radial profile of comet C/2007 N3 (Lulin) was plotted to assess the quality of the 
data for calculating a dust production rate. The radial profile of C/2007 N3 (Lulin) in the 



V band shows a deviation from the conical 1/p profile (IGehrz fc Ney 



19921 ). suggesting 



contamination from gas such as C2 at 5141 A. We therefore only calculate the dust 
production in / band. The / band radial profile of C/2007 N3 (Lulin) is shown in Fig. El 

To estimate the rate of dust production in comet C/2007 N3 (Lulin), we utilize the 



Afp quantity introduced by 



A'Hearn 



(119841 ) . This quantity serves as a proxy for dust 
production and when the cometary coma is in steady state, the value for Afp is an aperture 
independent parameter, 
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4 ^2 ^2 ^Q-O.4(m -m comet 

Afp= —± 



cm 



(1) 



where A is the Bond Albedo, / is the filling factor of the coma, m & is the apparent solar 
magnitude, m comet is the measured cometary magnitude, p is the linear radius of the 
aperture at the comet's position (cm) and and A are the heliocentric and geocentric 
distances measured in AU and cm, respectively. Cometary magnitudes are observed to 
follow similar phase angle effects as asteroids; therefore we also apply the phase angle 
correction of m comet (a = 0) = m comet (a) — Ca where a is the phase angle in degrees and 
C is correc tion factor o f 0.03 magnitudes per degree, the mean of the correction factors 



derived by 



Meechetal 



(119871 ). Figure H] illustrates the progression of Afp as a function of 



p and Table [2] reports values of Afp at a selection of distances from the comet photocenter. 



3.2. NIR and Polarimetry 



Comet C/2007 N3 (Lulin) was sufficiently distant at the epoch of our observations 
that we did not have the spatial resolution to map changes in polarization across the coma. 
Thus it is unknown whether or not there was any significant variation in polarization 
with coma morphology. Comet C/1995 01 (Hale-Bopp), for example, had a high surface 
brightness jet th at showed distinctly g reater fractional polarization in the NIR than the 



rest of the coma ( jJones fc Gehrz 



20001 ) . The only other comet with good spatial resolution 



and polarimetry at NI R wavelengths is comet 73P/Schwassmann-Wachmann 3 (SW-3), an 



ecliptic-family comet ( IJones et al. 



20081 ). SW-3 passed sufficiently close to the Earth to 



allow us t o measure the sur face brig htness and p o 



of 40 km (IJones et al. 



20081 ). While 



Jones et al 



arization of the coma with a resolution 



( 120081 ) and 



Harker et al. 



( 120111 ) found 



strong evidence for significant breakup of the dust aggregates released from the nucleus 
of SW-3 taking place over distances of 40-400 km from the nucleus, only small changes in 
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the fractional polarization across the inner coma were observed. All of the polarization 
observations of SW-3 were at phase angles > 20° and showed weak red polarimetric color. 
Despite the strong evidence for significant breakup of dust aggregates, the polarization of 
SW-3 is not distinguishable from other comets. 

The surface brightness profile of C/2007 N3 (Lulin) at H(1.65 /xm) is shown in Figure 
|5] along with the profile of a star for comparison. Given the seeing of 1'.'4 full-width 
half- maximum (FWHM), the surface brightness of C/2007 N3 (Lulin) is en tirely consisten t 



with a simple 1/p dependence as expected for a constant velocity outflow (IGehrz fc Ney 



19921 ) . The polarimetry was not of sufficient signal-noise to determine if there were any 



variations in fractional polarization at large distances (> 1000 km or 3") from the center. 



There is very little IN 



C/1995 01 (Hale-Bopp) flJones fc Gehrz 



IR polarimetry of co mets at phas e angles < 20°, pr imarily 



single observation of 65P/Gunn (jKelley et al 



200 01), lP /Halley (Brooke et al 



19871 ). and a 



20051 ) . A compilation of NIR polarimetry of 



comets is shown in Fig. El The solid line in this figure is the mea n phas e dependence of 



the polarization in the R band found by 



Levasseur-Regourd et al. 



(119961 ) . The majority, 



but not all, of the NIR polarimetry lies above this line, at least out to phase angles of 60°, 
consistent with the generally red polarimetric color of comets. The precise (low formal 
error) observat ions of comet C/1995 01 (Hale-Bopp), which contained unusually small 



dust particles (jWilliams et al. 



1997 



Mason et al 



2001 



Harker et al 



2002 



Wooden et al. 



19991 ). shows no negative branch. The lack of a negative branch in the NIR for C/1995 01 
(Hale-Bopp) can, to first order, be explained by the added polarization of light scattered 
from more 'Rayleigh-like' particles coming from the jet. The addition of these smaller 
particles increases the fractional polarization above what would otherwise be a typical 
polarization vs. phase curve for the rest of the coma. 



Our H band polarimetry of comet C/2007 N3 (Lulin) is the first observation of the 
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polarization of a comet in the NIR at a phase angle within a few degrees of zero. The 
polarimetry for comet C/2007 N3 (Lulin) given in Tabled] is plotted vs. phase angle in 
Fig. [71 The solid line is the typical dependence of polarization on phase angle in the R 
band mentioned earlier. Our R band observations of comet C/2007 N3 (Lulin) are entirely 
consistent with the typical optical behavior of polarization comets. The NIR polarization 
of comet C/2007 N3 (Lulin) is also entirely consistent with the contemporaneous visual 
polarization measurements, and the NIR negative polarization is clearly comparable with 
a typical optical negative branch polarization behavior observed in most comets at small 
(< 25°) phase angles. 



3.3. The 10 fim Silicate Feature 



Within the 10 /zm spectral range covered by the Spitzer SL1 module, the observed 
SED of comets often are comprised of amorphous carbon grains, which produce the 
underlying featureless emission (continuum) in the 8 — 13 fim wavelength region, and small 
( < 1 um) siliceous dust grains which prod uce broad features and d i stinct resonances in 



exces s of the continuum (for a reviews see 



Hanner Zolensky 



2010 : 



Kelley fc Wooden 



20091 ). Amorphous silicates with chemical composition (stiochiometry) similar to olivine 
(Mgj / ,Fe M ) 2 Si0 4 and pyroxene (Mg I ,Fe (1 _ x . ) )Si0 3 with x = y = 0.5 (Mg/(Mg+Fe) = 0.5) 
reproduce the broad width of the 10 /im feature. The distinct 9.3 and 10.5 Aim emission 



features are attributed to Mg-rich orthopyroxene ( jWooden et al. 



1999 



Harker et al. 



20021 ). 



while Mg-rich crystalline olivine is uniquely i dentified through its distinct, relatively narrow 



11.2 fim silicate feature (IHanner et al. 



grains with 0.9 < x ~ y<1.0 ( [Wooden 



1994 ) . Mg-rich crystallin e 



2008 



Koike et al. 



2003 



species are de 


ined 


Chihara et al. 


2002) 



Frequently in comets, the Mg content of crystalline silicates is often significantly larger than 
that of the amorphous silicate grains. 
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Weak 10 /im silicate emission is evident in the IRS spectra of comet C/2007 N3 (Lulin) 
(Fig. |H]). At 10.5 /im the silicate emission is 8.4 ± 0.1% above a blackbody curve fit to 
continuum points around 8 /im and 12.5 /im. The best fit blackbody has a temperature 
of ~ 228 K. The silicate feature strength at 10.5 /im of C/2007 N3 (Lulin) is relatively 



weak com pared to other comets, most ecliptic-fami. 



strengths (IKelley fc Wooden 



2009 



SitkoetaL 



200J). 



y comets have > 15% silicate feature 



he 10 /im SED of comet C/2007 N3 (Lulin) was modeled using the lHarker et al. 
(120021 ) thermal emission dust code which assumes that a collection of optically-thin, discrete 
(singular mineralogy) dust particles reside at the heliocentric (r^) and geocentric (A) 
distance of the comet at the ep och of observations (2008 Oct 04.3 UT) and adopts a Hanner 



grain-size distribution (HSGD; 



Hanner et al. 



19831 ) for n(a) da. The 10 /im mineralogy 



used in the model is derived from labo ratory studies of in terplanetary dust particles 



(Wooden et al 



(IBrownlee et al 



2000), micrometeorites ( IBradley et al. 



thermal models (e.g.. 



19991 ) . the NASA Stardust mission 



2006) as well as other 



Hanner et al. 



1994 : 



rain species employe d in other remote sensing 



Wooden et al. 



20041 ) . Additional details of the our 



therm al modeling specifically of comet C/2007 N3 (Lulin) is discussed in 
(2011). 



Woodward et al. 



The best-fit thermal models (Fig. [HI and given in Table [3]) suggest that the silicates 
in C/2007 N3 (Lulin) are dominated by pyroxene grains, and include minor amounts 
of crystalline olivine and orthopyroxene. Grains in the coma of C/2007 N3 (Lulin) are 
relatively large, the HGSD peaking at a p = 0.9 /im, and moderately porous (fractal 
porosity parameter D = 2.73) with a large grain slope N = 4.2. The submicron sized 
silicate-to-carbon ratio inferred from our models is 0.48 ± 0.06. The error bars on 
parameters derived from the thermal modeling reflect a 95% confidence limit. We also 
find that the silicate crystalline mass fraction for the submicron to micron-size portion 



-14- 



of the grain size distribution (IHarker et all 120021 : iMoreno et al.ll2003l ). defined as f cr yst = 
(crystalline)/(crystalline + amorphous), for comet C/2007 N3 (Lulin) is 0.14 ± 0.04. In 
contrast, other NIC comets such as C/1995 01 (Hale-Bopp) which had a very strong 
10 /im silicate feature had coma dust gr ains that were fract ally very porous (D = 2.5) 



with the HGSD peaking at a p = 0.2 /im (IHarker et al 



20021 ) and mineralogically diverse, 



i ncluding amorphou s and crystalline forms of both olivine and pyroxene with /, 



( IHarker et al 



crys 



0.68 



20021 ). The dynamically new comet C/2001 Q4 (NEAT) which exhibited a 
modest 10 /im silicate emission had dust grains that were 'solid' (D = 3.0) with a HGSD 
peaking at a p = 0.3 /im and a with f crys t = 0.71 ( IWooden et all 120041 ) . 



In Fig. [9] we plot the strength of the 10 /im silicate emission above the continuum, 
F w /F c , where F w is integrated feature flux over a bandwidth from 10 — 11 /im and F c is 
that of the local continuum at 10.5 /im, vs. the color temperature excess above black body 



equili 

(e.g., 



jrium at the dist ance of C/2007 N3 (Lulin), and various other comets, from the Sun 



Sitkoetal 



20041 ). The color temperature excess is defined as: 



T, 



Fit 



T, 



BB 



T, 



Fit 



278/Vr h 



(2) 



where Tp it is the black body fit to the observed continuum and is in AU. For comet 
C/2007 N3 (Lulin) this temperature excess is 1.133 ± 0.051, while F l0 /F c = 1.084 ± 0.011. 
Comet C/1995 01 (Hale-Bopp), which had a jet with much smaller grain aggregate or 



200(1 1. 



monomer size dust particles and no negative branch in polarization ( jJones &: Gehr: 
shows an exceptionally strong silicate feature and a large color temperature excess, also 
i ndicative of small gr ain aggregates or small grain monomers separated from the aggregate 



(I Williams et al. 



19971 ). C/2007 N3 (Lulin), in contrast, has a very weak silicate feature 
with little color temperature excess. The large compact aggregates that are required by 
the polarimetry are consistent with the thermal IR measurements that show a very weak 
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silicate feature. 



The silicate feature strength and continuum temperature e xcess of fragments B and 



C of co met 73P/Schwassm ann-Wachmann observed by Gemini ( iHarker et al.l 120111 ) and 



Spitzer f lSitko et al. 



201 ll ) are also included in Fig. [9] for completeness. For fragment B, the 



silicate feature strengths and temperature excesses were 1.213 ± 0.015 and 1.066 ± 0.014 
(Gemini), and 1.242 ± 0.002 and 1.177 ± 0.001 (Spitzer) respectively. For fragment C, the 
silicate feature strengths and temperature excesses were 1.314 ± 0.020 and 1.125 ± 0.020 
(Gemini), and 1.335 ± 0.003 and 1.158 ± 0.002 (Spitzer) respectively. 



3.4. Aggregate Models and Comet NIR Polarimetry 



The negative polarization branch, where the maximum linear polarization is of order 



2.0% is evident in all comets observed at visual waveleng 



(for a review see 



Mishcnenko et al. 



2010 



Kolokolova et al. 



hs at small, ^$25°, phase angles 



20041 ) . As we mentioned above, 



computer simulations of this characteristic, based on aggregate models, requires either large 
aggregates (thousands of monomers) or monomers made of rather transparent material 
(silicates). These simulations are usually based on ballistic aggregation that not only 
allows fitting the observational data but are a l so consisten t with the origin and evolution 



of comet ary dust. For example, 



Kimura et al 



( 12006 



20031 ) considered two types of ballistic 



aggregates: BPCA (Ballistic Particle Cluster Aggregates) and BCCA (Ballistic Cluster 
Cluster Aggregate). These two types of aggregates differ in their porosity. For an equivalent 
number N of constituent monomers of radius a, the BPCA aggregates are more compact 
than BCCA, where the porosity is defined as P = 1 — N(a/a c ) 1 ^ 3 and t he characterist ic 



radiu s a c = a/5/3 a g , where a g is the gyration radius of the aggregate (jKozasa et al 
1992 ). a g = ( 1/2N 2 ) x ^f^J r ji) — r(j)) 2 , with r(i) th e posi tion of the center of the 



:th 



monomer. 



Kolokolova et al 



(120071 ) and 



Kimura et al. 



(120061 ) describe the polarimetric 



- 16 - 



and IR properties of comet dust using these two types of ballistic aggregates, consisting 
of sub- micron monomers (of radius a ~ 0.1 fim) with a Halley-type composition which 
includes silicates, amorphous carbon, and organic refractory material. These models can 
account for the general behavior of the maximum polarization, the shape of the polarization 
curve as a function of phase angle P(a), and the negative polarization branch. In addition, 
such models also yield the low geometric albedos (of order ~ 5% t hat is typical for come t 



dust) and red photometric colors (normalized reflectivity gradient; 



Jewitt fc Meech 



that are typical for comet dust, for examp 



others discussed in 
Indeed, Fig. 2 of 



Kolokolova et al 



(2004 



e C /2004 Q2 (Machholz) ( 



or 



Lin et al 



Hadamcik Levasseur-Regourdl (120091 ) 



1986|) 



2007) and 



Joshi. Ganesh. fc Baliyanl ( 120111 ) indicates that the colors of comet 



C/2007 N3 (Lulin) are not blue. For proper characterization of the observed low geometric 
albedos and red photometric of com ets, the values of the refractive index used to describe 



the aggregate ensembles are crucial. 



Kimura et al. 



(120061 ) demonstrate that if the imaginary 



part, k, of the index of refraction becomes smaller than 0.4 the color of comet dust becomes 
blue. 



The same ballistic aggregates, BPCA and BCCA, were considered by 



Lasue et al. 



( 20091 ) to model comet polarization. Mixing them with spher oidal silicate 



Lasue et al. 



partic les they 



(I2009h do not 



achieved a good fit to the polarization phase curve. However 
provide any photometric characteristics of the dust, specifically the albedo and photometric 
color. This omission does not permit a critical examination of their model t o assess its 



validit y in the study of comet C/2007 N3 (Lulin). The comprehensive study of 



the 



Kimura et al. 



(120061) indicates that o ne might expect, due to the rather low absorption properties of 



Lasue et al. 



(120091 ) materials (they considered the r efractive index of o rganics equal 



to 1.88 + iO.l and of silicates 1.62 + % 0.003), that the 



Lasue et al. 



((2009J) dust would 



demonstrate rather high albedo and blue c olor. The erroneou sly high albedo and blue 



dust color is also apparent in the results of 



Shen et al. 



(120091 ) which uses aggregates to 
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model comet dust; a aggregate model first developed by IShen et all ( 120081 ). In the latter 
work, they consider random ballistic aggregates (BAs) constructing BAM1 and BAM2 
particles describing differing monomer migration after randomized first contact aggregation. 
BAs have porosities similar to BPCAs, yet are le ss porous than BCCAs . However, the 



2009 



20081 1 models are low, 



y which lead to r ather 



complex refractive indices of the silicates used in ( IShen et al. 
1.72 + i0.029 and 1.71 + i 0.031 in the visible and IR respectivel 
high albedo and erroneous photometric colors. For instance, the 

produces an albedo of order ~ 12% and blue dust color. Currently, the model m ost 
consistent with the observational and in situ data is that of iKolokolova fc Kimural (120 lOf ) 



Shen et al 



(120091) model 



where, following the findings o f Stardust mission, a m i xture of aggregated and solid particles 



was considered. The results of 



Kolokolova fc Kimural (120101 ) modeling are characterized not 



only by the correct polarization and photometric properties of the dust, including albedo 
and colors, but also by the model organics-to-silicate ratio and ratio of aggregates to solid 
particles that are consistent with the results of in situ studies of comets. 

An important constraint in interpreting polarimetric aggregate models applicable to 
our study of C/2007 N3 (Lulin) near perihelion is our 10 /im spectra^ which show the 
comet near the epoch of the optical/NIR polarimetric observations had a weak silicate 
feature (§ 13. 3p . Thermal modeling of the 10 /im SED suggests that the peak of coma grain 
size distribution a p is of order 1 /im with a large grain slope indicative of a population 
tail of grains are extant whose sizes are >> 1.0 /im. These lat ter grains are large as well 
as moderately porous, Pthermai ~ f vacuum = 1 - (<V<W a o) D ~ 3 ( Marker et all 120021 ) where 



a Q = 0.1 /im and D = 2.73. Large aggregates are also necessary to explain the negative 
polarization in the NIR and the red color of light scattered by comet C/2007 N3 (Lulin) 's 



3 The synthetic apertures used to determine the polarization (see Table d]) are comparable 
to the Spitzer extraction apertures. 
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dust. Extensive optical photometry of the comet near opposition (IJoshi. Ganesh. fc Baliyan 



20 111 ) also indicate a dominance of grains in the coma larger than 0.1 /im. Thus, the results 
derived from the scattered light and the thermal emission spectra are consistent, and the 
porosity of grains from our thermal models ( > 70% for a gra i n = 10 /im) suggests that the 
coma of comet C/2007 N3 (Lulin) is composed of grains with properties more like BPCA 
grains rather than BCCA grains. 

Negative polarization was not observed in the NIR in comet C/1995 01 (Hale-Bopp) 
whose dust was characterized by small particles. Disappearance of the negative polarization 
in this instance signals that at wavelengths > 1.0 /on, the size parameter of the dust particles 
became smaller than the wavelength (Rayleigh particles) for which polarization is always 
positive. This is not the case observed in comet C/20 07 (N3) Lulin. The genesis of negativ e 



Petrova. Tishkovets. &: Jockers 



20071) 



polarization is complex (for a detailed discussion see 
It results from a combination of properties of individual monomers, multiple scattering, 
coherent backscattering and near-field effects. At some monomer sizes these latter effects, 
combined with porosity and refractive indices, tend to produce effects that work in the same 
direction increasing the negative polarization, but at other monomer sizes they produce 
opposite effects. This is why one can get equal negative polarization for different types of 
aggregates or a big change in ne gative polarization for a sma ll change of, for example, the 



size parameter of the monomer (IKolokolova Kimura 



20101 ) . Computation of negative 



polarization is still a modeling challenge deserving future attention. 

In our analysis of comet C/2007 N3 (Lulin), we claim that the negative polarization 
does not contradict conclusions regarding the particle properties inferred the red polarization 
color (from the i?-band through the iJ-band) and those derived from the thermal IR 
observations that indicate the grains are large (i.e., micron-sized or larger) and of low 
porosity (perhaps as compact as an aggregate of submicron particles can be). Another 
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explanation for the negative polarization observed in the coma of comet C/2007 N3 (Lulin) 
is that the dust was composed of transparent silicate particles with little or no carbon 
content. However, cometary dust is usually characterized by a low albedo that presumes 
a large content of carbon, and photometrically C/2007 N3 (Lulin) does not show any 
peculiarity that would allow us to suppose that its albedo is significantly different. Also, in 
the case of abu ndant silicate part icles, we would expect a noticeable spectral dependence 



of polarization (IZubko et al. 



20091 ). In addition, our thermal model has a silicate to carbon 



ratio of 0.48, inconsistent with dust dominated by transparent silicates. 



CONCLUSIONS 



We have found that C/2007 N3 (Lulin) clearly has a typical optical negative branch 
in the polarization as well as exhibiting a negative branch in the near-infrared at low 
phase angles. When compact aggregates are larger than the longest wavelength considered 
it is known that the depth of the negative branch does not depend significantly on 
the wavelength; furthermo re, the wavelength d ependence becomes less pronounced as 



the phase angle decreases (IBelskaya et al. 



20091 ). We have also found that these large, 



low-porosity aggregates in comet C/2007 N3 (Lulin) are consistent with the thermal infrared 
measurements that showed a very weak silicate feature. 

Our thermal model contains moderately porous grains (D = 2.727), with a peak 
grain size of 0.9 /zm. This is qualitatively consistent with the large aggregate grains 
needed to explain the polarimetry if the monomers are ~ 0.1 /im in size. We conclude 
that the dust in comet C/2007 N3 (Lulin) is dominated by large and compact aggregate 



particles, made u p of thousands of smal 



periodic comets (IKolokolova et al. 



porous aggregates. However, 



monomers. Compact aggregates are typical for old 



2007) but no t for new ones, which tend to have more 



Sitko et al 



(120041 ) show that the dust properties of some 



-20 - 



nearly-isotropic (including Oort cloud comets) and ecliptic comets are similar based on 
intercomparison of the strength of the 10 /j,m silicate feature versus the color temperature of 
the dust. Our polarimetry of comet C/2007 N3 (Lulin), an Oort cloud comet with thermal 
emission properties equivalent to ecliptic-famil y comets, demon strate that grain structure 



may account for this observation (see Fig. 4 of 



Sitko et al 



20041). 



The occurrence of a negative branch in the polarization in the near-infrared is likely 
typical for most comets. Comet C/1995 01 (Hale-Bopp), which showed no negative branch 
in the near-infrared, must be considered an anomalous case due to the production of 
significant numbers of very small, submicron grains. These latter small coma grains could 
be individual monomers that have broken off from the larger, porous aggregates. For 
C/2007 N3 (Lulin), the more compact aggregates must maintain most of their integrity 
after release from the nucleus. Although aggregate models for the polarization and the large 
porous grain model for the thermal emission of comet C/2007 N3 (Lulin) share qualitative 
features, a rigorous quantitative gr ain/grain-aggregate model t hat is self consistent for 



both scattering and emission (e.g., 



Kolokolova fc Kimura 



2010) has yet to be applied to 



observations of individual comets. Such an effort will be an important future work. 
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Table 1. POLARIMETRY OF COMET C/2007 N3 (Lulin) 



Date (2009) 


UT (Hr) 


Phase (deg) 


Filter 


P(%) 


21 Feb. 


8:00 


21.3 


R 


0.41 ±0.15 


22 Feb. 


7:30 


17.4 


R 


-0.50 ±0.15 


26 Feb. 


10:00 


0.44 


R 


-0.40 ±0.15 


27 Feb. 


7:30 


4.0 


R 


-1.30 ±0.15 


27 Feb. 


7:30 


4.0 


H 


-1.20 ±0.15 



a i?-band values determined using a 12" diameter synthetic circu- 
lar aperture, while the if -band value is measured in a 7'.'5 diameter 
synthetic circular aperture. 



Table 2. Af p VALUES FOR COMET C/2007 N3 (Lulin) 



Aperture [arcsec] 


p [km] 


I [mag] 


Af p [cm] 


28.35 


10000 


8.49 ± 0.06 


3188 ± 111 


42.3 


15000 


8.05 ± 0.07 


3213 ±138 


70.2 


25000 


7.49 ± 0.07 


3228 ± 138 


281.25 


100000 


6.24 ± 0.07 


2545 ± 109 
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Table 3. BEST-FIT THERMAL MODEL AND DERIVED PARAMETERS 



Sub-/iiii 



Dust component 


N p x 10 16 b 


mass fraction 


Amorphous pyroxene 
Amorphous olivine 
Amorphous carbon 
Crystalline olivine 
Crystalline pyroxene 


1092^6 

5814+^ 
2681™ 
2481^ 


0.17 ±0.01 
0.015 ±0.005 
0.68 ±0.03 
0.07 ±0.02 
0.07 ±0.04 


Other model parameters 


Degrees of freedom 
Total coma mass 
Silicate / carbon 


47.2 
101 

(1.096 ± 0.048) x 10 5 kg 
0.48 ±0.06 




Uncertainties represent the 95% confidence level. 



b Number of grains at the peak of the grain size distribution. 
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Fig. 1. — Positions of the Earth (triangles) or Spitzer (upside-down triangle), and comet 
C/2007 N3 (Lulin) (filled circles) on: (1) 2008 October 04 UT, (2) 2009 February 27 UT, 
and (3) 2009 March 03 UT depicted in heliocentric ecliptic coordinates. The orbits of the 
terrestrial planets, Mercury through Mars, and comet C/2007 N3 (Lulin) are drawn as gray 



lines, and the position of the sun indicated by the 0. The axis tick marks are in units of 
AU. 



Fig. 2. — (top) Greyscale V band (left) and / band (right) images of comet C/2007 N3 
(Lulin) from the Bok 2.3-m Telescope, (bottom) same as above with isophotes overlain on 
the image from 17 to 19 magnitudes per square- arcsecond at 0.25 magnitude per square- 
arcsecond intervals. The field of view in each filter is 18.4 arcminutes x 14.3 arcminutes. 
The grey scale intensity display table is proportional to the square of the counts. 
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10 2 10 3 1 4 10 5 10 6 

Rho [km] 



Fig. 3. — Radial profile of fluxes as a function of linear radius as measured in / band from the 
photocenter of comet C/2007 N3 (Lulin) obtained on 2009 March 03 UT. Solid line denotes 
a 1/p profile. 
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Fig. 4. — Afp as a function of p derived from / band observations of comet C/2007 N3 
(Lulin) on 2009 March 03 UT. 
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Fig. 5.— An E-W cut across comet C/2007 N3 (Lulin) in the H (1.65 /jm) band on 2009 
February 27 UT. The profile of a star is shown as a dotted line. The solid line is a 1/p fit to 
the comet coma profile. 
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6. — Compilation of comet polarimetry in the near-infrared (adopted from 



Kelley et al. 



2004 f) . The solid line is the a verage trend for comets in the visual R (0.676 /im) band 



( iLevasseur-Regourd et al. 



1996J). The H (1.65 /im) measurement for C/2007 N3 (Lulin) 



from this work is inset as the blue symbol. 
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Fig. 7. — Polarization vs. phase angle for C/2007 N3 (Lulin) from data in Table [TJ 
he solid line is the average behavior for comets in the visual R (0.676 /mi) band 



( ILevasseur-Regourd et al 



19961 ). The polarimetry of C/2007 N3 (Lulin) is entirely consistent 



with this trend both at R and at if (1.65 jttm). 
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Fig. 8. — Spitzer IRS 7-13 fim spectrum of comet C/2007 N3 (Lulin) obtained on 2008 
October 04.3 UT (r^ = 1.90 AU; Spitzer-comet distance = 1.67 AU) derived from an 
extraction aperture of 9.25" x 9.25" centered on the nucleus. Data are shown in black dots, 
(top) The composite thermal model best-fit spectral energy distribution (in AFa vs. A space) 
is indicated by the solid red line, while the contribution from amorphous pyroxene (blue), 
amorphous carbon (orange), and crystalline olivine (green) are indicated, (bottom) The 
observed IRS flux density divided by a ~ 228 K blackbody continuum (Fa/Fa, t) to highlight 
spectral details of the 10 /im silicate feature. The 'feature' near 7.5 /zm is an artifact of IRS 
order overlap. 
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Fig. 9. — The strength of the 10 /xm silicate emission feature, F w /F c , as a function of the 
excess color temperature above the equivalent blackbody equilibrium temperature at th e 
Sun-comet distance, T c /Tbb- The grey points are all comets listed by 



Sitko et al. 



( 120041 ) . 



The individually labeled are c omets C/2007 N3 



Wachmann fragments B and C (IHarker et a 



01 (Hale-Bopp) observed at = 2.73 AU (ISitko et al. 



2011 



ulin) from this work, 73P/Schwassman- 



Sitko et al. 



20111), lP/Halley, and C/1995 



2004J). These 5 latter comets all have 



measured mic - 



ni 



Sitko et al. 



R sp ectra and near-IR polarization. The trend-line is the same as presented 



(|2004j) 



